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Abstract 

This article is the continued version of the analytical solutions for the pressureless Navier- 
Stokes equations with density-dependent viscosity [9]. We are able to extend the similar 
solutions structure to the case with pressure under some restriction to the constants 7 and 9. 
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1 Introduction 

The Navier-Stokes equations can be formulated in the following form: 

pt+V ■ (pu) =0, 

(1) 

(pit) t +V • (pu (g) u) + 5X7P =vis(p, u). 

As usual, p = p(x, t) and u(x, t) are the density, the velocity respectively. P — P(p) is the pressure. 
We use a 7-law on the pressure, i.e. 

P[p) = Kp\ (2) 

with K > 0, which is a universal hypothesis. The constant 7 = cp/c v > 1, where c p and c v are the 
specific heats per unit mass under constant pressure and constant volume respectively, is the ratio 
of the specific heats. 7 is the adiabatic exponent in ((2|). In particular, the fluid is called isothermal 
if 7 = 1. It can be used for constructing models with non-degenerate isothermal fluid, d can be 
the constant or 1. When S — 0, we call the system is pressureless; when 8 = 1, we call that it is 
with pressure. And vis(p,u) is the viscosity function. When vis(p,u) — 0, the system ([1]) becomes 
the Euler equations. For the detailed study of the Euler and Navier-Stokes equations, see p] and 
[3]. Here we consider the density-dependent viscosity function as follows: 

vis(p,u) = vG"(p) V - u )- 

where p(p) is a density-dependent viscosity function, which is usually written as p(p) = up 6 with 
the constants k, 8 > 0. For the study of this kind of the above system, the readers may refer [6], 
8, 

The Navier-Stokes equations with density-dependent viscosity in radial symmetry can be ex- 
pressed by: 

p t + up r + pu r + ^-pu = 0, 

(3) 

p (u t + uu r ) + VKpi = { K p e ) r (£=i u + u r ) + Kp e (u rr + ^u r - Z^u), 

Recently, Yuen's results [S] showed that there exists a family of the analytical solutions for the 
pressureless Navier-Stokes equations with density-dependent viscosity: 
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for 9 = 1: 



a W = a (°) = «o > 0, a(0) = ai, 

y( X ) = 2Nn x2 + Q > 



(4) 



for 6^1: 



0, fory(^ y )<0 
a(t) = a ^~Ne-N+2 , o(0) = a > 0, a(0) = oi, 



(5) 



= "-^(0-1)^2 + ^-1, 
where a > 0. 

In this article, we extend the results form the study of the analytical solutions in the Navier- 
Stokes equations without pressure [9] to the case with pressure. The techniques of separation 
method of self-similar solutions were also found to treat other similar systems in [2] , [I] , [5] , [7] , 
[S] and 0. 

Our main result is the following theorem: 

Theorem 1 For the N -dimensional Navier-Stokes equations in radial symmetry f3|), there exists 
a family of solutions, those are: 
for8 = j= 1, 

P{t,r) 



Ae Bi ^ft) ) + ° (. \ a(t) 



a(ty 



*(*) " + = °. o(0) = fl o > 0, a(0) = 0ll 



(G) 



where A > 0, B and C are constants, 
for 9 = 7 > 1, 



_ff 7 





SpF ^ a(t)««+i a(t) 



> u ^ r > = W) r ' 

0,fory(^)<0 
-^tw$s = 0, a(0) = a > 0, a(0) = oi, 



(7) 



y(z) = •-fa(e-l)z* + afl-i, 
where ao, a\ and a > are constants] 



4 



Yeung and Yuen 



J UI 2 ~ 2 /V /V 



o ; /«ry(^j)<0 

a(i) = <r(mi + n) s , < s = 7jV _ 2 jV+2 < 1 

^rry(^)^ 2 - ^y(^) e ' 2 ] y(z) = ( -^£z, y(0) = a > o 

where m, n > 0, cr > and a are constants. 



(8) 



2 Separation Method of Self-Similar Solutions 

Before we present the proof of Theorem [1] Lemmas 3 and 12 of [9] are quoted here. 

Lemma 2 (Lemma 3 of |9j) For the equation of conservation of mass in radial symmetry: 

N - 1 

p t + Up r + pu r H pu = 0, 

r 



there exist solutions, 

p(t r) = /(r/a(t)) „(tr)A 
Pl ' j o(t)* ' 1 ' ' a(t) ' 

wifft i/ie /orm / > S C 1 and a(t) > £ C 1 . 

Lemma 3 (Lemma 12 of |9j) For the ordinary differential equation 

y(z)y(z) n -& = 0, 
j/(0) = a>0,n^ -1, 
where £ and n are constants, 
we have the solution 



y( z ) 



* +1 'i(n + l)£z 2 + a™+ 1 , 



(9) 



(10) 



(11) 



(12) 



where the constant a > 0. 



At this stage, we can show the proof of Theorem [T] 
Proof. Our solutions ©, Q and |(5J| fit the mass equation ©i by Lemma ©. Next, for the 
equation ((3]) 2, we plug our solutions to check that. 
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For 9 = 7 = 1, we get 



:n-i \ . N-i 

p(U t + U ■ U r ) + K [p\ r - (Kpjr U + U r - Kp r {U rr H U r - 



--3 + K ^ B 

a{t) 6 a{t) 

Ae B(wj) 2 +C 



a {ty 



-2r 



ANne B( ^ )2+C 



-B 



a (ty 



a(t) 



2BK BNna{t) 
~ajf) + a{tj 2 



= 0, 



where the function a(t) is required to be 



, s 2BK BNn'ait) 
d(t) -— + j-^- = 0, a(0) = a > 0, a(0) = a%, 



a(t) a(t) 2 



where a and a\ are constants. 
For 6 = 7 > 1, we have: 



p{u t + u ■ u r ) + K [p J ] r — (np ti ) r ( u + u r ) — np e r (u rr H u r 



a(t) N a(t) a(t)'r N + 1 
By defining z := r/a(t), and requiring 



a(t) 



a(t). 



y(z)z~y(z) e ~ 1 y(z) = 0, 
z-y(z) e - 2 y(z)=0, 



T9l) becomes 



y(z)z 



a(t) K9 Nn6a{t) 



a(t) N a(t) 0N + 1 a(tY 
where the function a(t) is required to be 

a(t) K~i Nn9a(t) 

+ a(*) 9Ar+1 ~ a(t) eN+2 

Therefore, the equation ^ 2 is satisfied. 

With n := 8 — 2 and £ := 1, in Lemma [31 (|2"T]) can be solved by 



= 0, 



= 0. 



z/0) 



1 1 



l)z 2 +a e -\ 



() 
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where a > is a constant. 

For the case of ^ + 5 — jj = > 1 ~ 7^, we have, 

f N - 1 

p(w t + u • u r ) + K [p 7 ] r - (np 1 ) r u + u r I — Kp° r (u rr + 



r?" H t^tt-ttti — — — a(i). 



AT-1 AT - 1 



-u) (25) 



a(i) w o(i) o(t)T*r+i a(t) aN+2 

By letting a(i) = a(mt + n) s , we have 

r S ( S -l)(mt + ny^m 2 ar KOyj^- 1 ^) amg jV K fly(-fa)«-iy(-fa) 

o-^mf + n)'* o(t) o-7JV+i( m f + n ) S ( 7 JV+i) «r«^+2( mi + n )»(fljv+2) 



(26) 



(mi + 



(27) 



r s(s — l)m 2 



^a(t)^ a N - x a{t) (mt + n )«*-(«-2) 1 CT 7Ar+1 (m< + n)"(T JV + 1 ) c^+^m* + „)«(WV+2)-(«-i) 



(28) 



r . s(s — l)m 2 r 



a(i)' ct^-1o(*) (m* + n) sJV " s + 2 cr^+^mi + n) 5 ^^ 1 ) a^+^mt + n) eN + s + 1 



Here we require that 



sN-s + 2 = s(jN+l), 
s(jN + 1) = s(9N + 1) + 1. 



That is 



< s 



(7 - 0)iV 7^V - AT + 2 



< 1. 



(29) 



(30) 



(31) 



In the solutions ©, it fits to the conditions (pTTj) by setting f + 5 — jt ~ > 1 — > and 



7JV-JV+2 



Additionally by defining 2 := r/a(t), the equation (|2T)|) becomes 



(mt + n) Ws - ;5 + 2 
Here we require that 

at 7 



(s — l)m 2 AT 7 



n-— z + ^n+M 2 ) y( z )-zeNTTy^ 



(32) 



7-2 



tuNkO 



_ S(J jN+l 

The proof is completed. ■ 

In the corollary can be followed immediately 



y{'- 



2/0) = 



(1 - s)m 2 



(33) 



Corollary 4 For m < ; the solutions Q) blowup at the finite time T = 



—m/n. 
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